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We show that hidden hot dark matter, hidden-sector dark matter with interactions that decouple 
when it is relativistic, is a viable dark matter candidate provided it has never been in thermal 
equilibrium with the particles of the standard model. This hidden hot dark matter may reheat to 
a lower temperature and number density than the visible Universe and thus account, simply with 
its thermal abundance, for all the dark matter in the Universe while evading the typical constraints 
on hot dark matter arising from structure formation. We find masses ranging from ~3 keV to 
~ 10 TeV. While never in equilibrium with the standard model, this class of models may have 
unique observational signatures in the matter power spectrum or via extra-weak interactions with 
standard model particles. 

PACS numbers: 95.35.+d, 98.80.Cq, 14.80.-j 



Introduction - Dark matter is a known unknown. While 
we know its abundance [l[ and some of its cosmological 
and astrophysical properties (see, e.g., Ref. Q) its funda- 
mental nature remains elusive. Viable candidates, that 
may also play a role in the solution of other pressing is- 
sues in physics, have been proposed including weakly in- 
teracting massive particles (WIMPs) 0], and axions 
Although they have drastically different properties both 
would have been relatively cold dark matter (CDM) par- 
ticles (relative to their mass) when first produced. They 
would have moved only a short distance, cosmologically 
speaking, before their initial momentum redshifted away. 

On the contrary, while standard model (SM) neutrinos 
have mass 0, they fail as a candidate for dark matter 
because they decoupled from the primordial plasma at 
a high temperature relative to their mass. They are hot 
dark matter (HDM) particles that stream through the 
Universe at velocities v ~ c and travel huge cosmologi- 
cal distances before slowing down. This free-streaming 
would erase small-scale density fluctuations in violent 
conflict with observations 0. Moreover, to account for 
a cosmological dark matter density of fldh 2 ~ 0.11 [l| 
a single Dirac (Majorana) fermion that decoupled while 
relativistic must have a mass ~ 50 eV (~ 100 eV) if it 
decoupled prior to the QCD and electroweak phase tran- 
sitions, and a lower mass if (like neutrinos) it decoupled 
later. As the Tremaine-Gunn bound on thephase space 
density of dark matter requires m > 300 eV [7| , it appears 
that an interacting species once in thermal equilibrium 
that decouples when relativistic, like the SM neutrinos, 
can only account for a minor fraction of the dark matter 
in the Universe. 

However, there is an important caveat to this conclu- 
sion. It assumes that dark matter was, at some point in 
the history of the Universe, in thermal equilibrium with 
standard model particles. This need not be the case if 
dark matter exists in a hidden sector (HS). In fact, dark- 
matter could have been thermalized with hidden sector 
particles at a HS temperature less than the SM tempera- 



ture. 1 This may occur straightforwardly, for instance, if 
reheating after inflation proceeds more efficiently to the 
SM sector than the HS. 2 Hidden sectors have long been 
imagined Q, and HS dark matter candidates that freeze 
out when non-relativistic have recently received attention 
Q. Potential implications of extra HS interactions have 
also been considered (see, e.g., Refs. [Ioj|). 

In this Letter we show that hidden hot dark matter 
(HHDM), a thermalized HS species with interactions that 
decouple while relativistic, is a candidate for the all of the 
dark matter in our Universe. This alternative class of 
dark matter candidates has a thermal relic abundance, 
significantly different particle properties than standard 
thermal relics, and might have extra-weak interactions 
with the SM. 

Mass and Relic Density - We assume that dark matter 
is a particle \ of mass m x which freezes out when it is 
relativistic. This analysis is for a fermion, but appropri- 
ately adapted results hold for bosons. A fermion with d x 
internal degrees of freedom in thermal equilibrium at a 
temperature T x has a number density 
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where T is the SM temperature, g s {T) is the effective 
number of entropy degrees of freedom coupled to the SM, 
and £ = T x /T is the dark-to-SM temperature ratio at the 
temperature T x . When HS interactions of x particles 
decoupled at temperature T x f we have £/ = T x f/Tf. 
Here £(3) ~ 1.202 and for Dirac fermion d x — 4 counting 
both particles and antiparticles (for a Majorona fermion 
d x = 2). If x particles are dark matter they are now 



1 A different possibility, one example being the Dodelson- Windrow 
sterile neutrino [Tl| , is non-thermal dark matter that never equi- 
librates with any species. This alternative is not our focus here. 

2 The relative temperature could also be reduced if g" decreases 
by a huge factor before electroweak but after hidden decoupling. 
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non-relativistic and have an energy density 
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where n x f is the \ number density at freeze-out, aj 
the scale factor at freeze-out and ao the scale factor at 
present. This means that x particles must have a mass 
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in order to be dark matter. Here g s ^ and <?q are the effec- 
tive number of entropy degrees of freedom in the Universe 
at x decoupling and at present respectively. 

We emphasize that as x particles decouple when rel- 
ativistic, the dark-matter density is determined by the 
thermal abundance of a relativistic species, and is ex- 
pressed in terms of the temperature and thermal his- 
tory of the Universe. Unlike non-relativistic freeze-out, 
it is independent of the interaction cross section except 
for requiring it be small enough to allow relativistic de- 
coupling. Instead we find, at fixed fldh 2 , there is a 
unique relationship between mass m x and the dark-to- 
SM temperature-ratio £/. 

As it is cold relative to the SM plasma, we show below 
the lower limit on the dark-matter mass from structure- 
formation are slightly relaxed and m x ~ 3 keV is vi- 
able. For small £/, this HDM might still be very heavy 
and, as discussed below, may have masses as large as 
m x ~ 10 TeV for the class of models we discuss. On 
the low mass end HHDM acts like warm dark matter 
and alters the dark-matter power spectrum on relatively 
large length scales where it is sensitive to constraints from 
cosmology or observations of substructure. Furthermore, 
while the thermal history of dark matter is determined 
by HS interactions this class of dark-matter candidates 
may also couple to the SM with extra-weak interactions. 
Extra-weak interactions can result, for instance, in the 
decay of dark matter to SM particles on timescales longer 
than the age of the Universe and imprint a characteristic 
energy scale 3keV < m x < 10 TeV on any SM decay 
products. 

Decoupling - We discuss here our physical assumptions 
and determine constraints on the decoupling of x from 
the HS plasma so that it is viable HHDM. Instead of 
focusing on a specific theory, we imagine that the anni- 
hilation of x particles is mediated by a set of particles 
with typical coupling A and high-energy mass scale tha- 
For instance, if x were similar to a SM neutrino then A 
would be a typical weak-scale coupling, and m\ would be 
close to the mass of the W boson. While not completely 
general this setup encompasses a wide variety of models 
and has a natural analogy with the weak interactions of 
SM neutrinos. 

If x particles annihilate via the exchange of a relatively 
massive particle we can estimate the annihilation cross 



section on kinematic and dimensional grounds. If the x 
particles have momentum p 3> m\ then the annihilation 
cross section goes like a a ~ A 4 /s where s is the Lorentz- 
invariant Mandelstam variable. However, for particles 
with momentum m x < p < Tn\, like for SM weak in- 
teractions, the propagator of the mediator is dominated 
by 771a and a a ~ G A s, where Ga = A 2 /m A is an ef- 
fective four-fermion coupling constant. For p < m x we 
expect a a ~ G\m x , but further kinematic suppression is 
possible depending on the structure of the effective inter- 
action. 

For hot production we are in the intermediate regime 
where a a ~ G A s and so, on dimensional grounds, 
the thermalized cross section (times velocity) goes like 
(o~ a v) ~ G\T 2 where T x is the x particle temperature. 
In viable models T x must be less than the SM temper- 
ature T. In terms of the temperature ratio £ the total 
thermalized annihilation rate is 



F x T x 



(4) 



where we have used n x cx T x for a relativistic species. 
We can write the total thermalized cross section as 
(a a v) = N^Gj^T 2 , where N c h is the effective number of 
channels involved in number-changing interactions. This 
factor accounts for the number of channels via which x 
can annihilate, differences in the true rates of the rele- 
vant processes, and factors arising from the thermal rate 
calculation. The thermalized coupling F x and effective 



coupling Ga are related by F x 



m^y/djiGA. 



The freeze-out temperature Tf is the standard model 
sector temperature below which 



r a \T xf =£ f Tf < H\ T} 



(5) 



and number-changing interactions become ineffective 
compared to expansion. In a model with HS degrees of 
freedom the Hubble rate during radiation domination is 
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where g(T) and gh(£,hT) are the effective number of de- 
grees of freedom in the SM sector and HS respectively, 
and £h = 2/,/T is the hidden-to-SM temperature ratio 
(equal to £ for T > Tf). In viable models we have 
£/ < 0.1 and thus <^ g (unless gh is very large). 
This means the SM sector dominates the expansion of 
the Universe and we recover the standard form 
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We use this approximation in what follows. 

Using Eqs. @ and we find a freeze-out temperature 



T f ~ 16.5 
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FIG. 1: The viable HHDM (m x ,F x ) parameter space for m A = 10 GeV (left), m A = 100 GeV (middle), and m A = 1000 GeV 
(right) is shown for Q.dh 2 = 0.11, d x = 4, iV c h = 1, g s Dr — <?£>/, and a standard model thermal history in the SM sector. The 
darker shading shows the allowed region if the free-streaming constraint and recouping constraint are imposed (the latter with 
£ r = 1000). The lighter shading shows how the allowed region expands if the phase space density and four-fermion constraints 
are used instead. The top scale, and the associated dotted contours, shows the dark-to-SM temperature ratio at freeze-out £/. 
For fixed models with lower F x , that freeze out prior to the QCD phase transition, require a larger m x to achieve a dark 
matter density Q.dh 2 = 0.11 compared to models that freeze out later. Viable corners of parameter space exist over the range 
10 -4 GeV < mA < 10 11 GeV , while the maximum \ mass allowed is m x ~ 10 TeV near tua — 150 TeV. 



We must ensure freeze-out occurs when \ particles are 
relativistic. Requiring that T x f = > 3m x leads to 

the constraint 

( F x V rn x N7/ 8< /c m y 1/3 / g/ y/4 
Vl0- 4 GeV" 2 /V2McV/ ~ \0A J \83J V ; 

on F x m x 6 , where C m = (Qdh 2 / d x )(g S f / g$). This is 
shown in Fig. ^ as the curve labeled Decoupling. 

We have assumed here an effective four-fermion model 
is a faithful approximation to interactions for T x < toa- 
This imposes consistency constraints on the model for a 
given high-energy scale m\. Requiring T x f < m\/5 we 
find the constraint 

( F x™T V ™ x y 1/3 > Z (^lY 1/3 (9iY^ 

\10- 7 GeV- 1/2 y^MeV/ ^2 v 0.75y V99/ 

(10) 

on F x m^ 2 m x 1 ' 3 . This is shown in Fig. [1] as the curve 
labeled Four-Fermion. If A < 1/3, so the dimensionless 
coupling A of the high-energy theory is perturbative at 
relevant energy scales, we find an upper limit 

iWS^f^-am^. (") 

This is shown in Fig.[T]as the curve labeled Perturbatvity. 
These consistency constraints might be relaxed or altered 
if a model is specified in detail. 



Recoupling - If the x number density is given by its ther- 
mal distribution function the HS interactions discussed 
above (at least minimaly without extra high-energy in- 
teractions) must be strong enough that x particles can 
thermalize in the early Universe after HS reheating. This 
places a potential recoupling constraint on the model. 

At very high momentum p 3> m\ we assume the 
number-changing cross section scales as o~ a ~ A 4 /s. This 
means for T > m\ the annihilation rate is approximately 

T Q ~ F x 2 m 4 T x . (12) 

We thus expect the ratio T a /H peaks near T x ~ m\. 
For these interactions to thermalize x particles we require 
T xr = £ x T r , the dark-matter temperature when 

T a \ Tr > H\ Tr (13) 

is first satisfied, to be greater than m\. In particular we 
require T xr > £ r m\, where £ r is a buffer parameter to 
allow for the time it may take the HS to thermalize from 
its initial distribution. In what follows we set £ r ~ 1000, 
and the limit on F x m^[ 2 m x 1 ^ is 

/ f x mfm- 1/3 \ > ( W,t \- 1,Z [e7 (9r\V* 
^lO-SGcV- 5 / 6 J ~ V 0.77 J V 103 V107J ' 

(14) 

where K s r f = {g s r / 'g})(g s D f I 'dhr)- This is sri0wn in Fig- HI 
as the curve labeled Recoupling, and applies if no other 
thermalization mechanisms are present in the HS. 
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Structure Formation - For low masses \ particles appear 
to us like warm rather than cold dark matter and ap- 
propriately adapted constraints from cosmological struc- 
ture formation apply. After last interacting x particles 
stream unimpeded through the Universe with v ~ c 
until their momentum (which redshifts as p oc a -1 ) 
is small compared m x . This occurs at a HS temper- 
ature T x nr — m x /3.15 and a SM-sector temperature 
T NR ~ (m x /(3.15C/))(5//ffN R ) 1/3 ■ We estimate the co- 
moving free-streaming length is (see, e.g., Ref. flil ]) 



^FS 



where 



82 
~9 
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/ 42.5kpc-0.137 \ /5nr 
V ^fs r ■ tt m h 2 J V So 




Afs r - 42.5 kpc 



(16) 



d x -0.11 



and we define, here and elsewhere, g = g s ■ (g s /g) 3 . In 
Eq. (fT6|) . Ai?s R is the distance the particle free-streams 
when ultrarelativistic, and Xfs includes the distance tra- 
versed when semi-relativistic prior to matter-radiation 
equality. As x particles stream isotropically features in 
the initial density field smaller than Xfs are erased. Lim- 
its from measurements of the matter power spectrum in- 
dicate Xfs ^ 110 kpc which corresponds, up to logarith- 
mic corrections, to the limit m x > 2.9 keV. This is shown 
in Figure. [1] as the curve labeled Free- Streaming. 

Another limit on m x arises from the Tremaine-Gunn 
bound [3]. The maximum of the coarse grained x phase 
space density (PSD) in a galaxy today with core radius r c 
and velocity dispersion a must be less then the maximum 
of the microscopic PSD (the Fermi-Dirac distribution) 
after interactions freeze out. This leads to the constraint 



> 
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4 km/s 



1/4 



50 pc 



1/2 



(17) 



if the x distribution is well approximated by an isother- 
mal sphere. For example, for the ultra-faint dwarf 
spheroidal galaxy Leo IV we have m x > 0.7 keV |13l. al- 
though rigorously this bound is slightly weaker [14j . This 
is shown in Figure. [1] as the curve labeled Phase Space. 



Extreme Models 



The maximum m x allowed, 



rn. 



5 TeV 



10 3 



0.11 



4/5 




is at rriA ~ 150 TeV, at the confluence of the decoupling, 
recoupling, and perturbativity constraints. If the four- 



fermion condition is required instead, at toa — 400 TeV, 
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(19) 



Here, we have g r j = (jrigr/gf) 1 ^ 3 , and for a SM thermal 
history g r j ~ g f ~ g f — 106.75. We find viable models 
require m x < 10 TeV. Although 3keV < m x < 10 TeV, 
viable (m x ,F x ) regions exist over the wide range of the 
high-energy scale 100 keV < m A < 10 8 TeV. 

Standard Model Interactions - A key requirement of this 
scenario is that the HS never thermalized with the SM 
sector. This limits the coupling between the HS and the 
SM. The HS might interact with the SM via several chan- 
nels, and each constraint is best handled individually for 
each particular channel and model. 

We illustrate here a case where four-fcrmion interac- 
tions with SM fermions dominate and scale like 07 ~ A^ /s 
above an interaction scale mi and as aj ~ G 2 s below it, 
where Gj = Xj/mj is an effective coupling for SM-HS 
interactions. We write the interaction rate as 



r, = 



i,jeSM 



n z ((a^ x v) + (a^ xx v)) = FfT 5 , (20) 



when T < mj. Here, the sum extends over all chan- 
nels that produce at least one x particle, connect the 
SM to the HS, and Fj ~ Gi is the thermalized cou- 
pling. As Th < T we limit the HS production rate from 
SM particles Tj = Tsm—>hs and not the smaller rate 
Fhs^sm- When T > mi the interaction rate scales 
as Tj ~ F 2 m\T . To determine the abundance Y X M = 
n x /s of x particles generated from SM plasma we inte- 
grate the appropriate Boltzmann equation for these pro- 
cesses to find y x SM = / * dir/F e «. Here Y eq = n eq /s is 
the abundance of a species in equilibrium with the SM 
plasma. The ratio Tr/H peaks near T ~ mi and most x 
particles generated via SM-HS interactions are produced 
near this maximum. For the rates above we find 



n v 4 



VkeV/ 
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(21) 



If we require n x M /n x < e SM we find the constraint 



tp 3/2 

rim j 



10- 10 TcV~ 1/2 / V 10keV 



< 2 6 ' ^ l9 ' ndk2 4 



0.01 1000 0.11 dy 



(22) 



on Fim 3 j 2 m\! 2 . Note we have neglected potential reso- 
nant effects near T ~ m; which may occur in some mod- 
els and would lead to even more stringent constraints. 
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Interactions of this type connecting the HS and SM 
must be extra weak, much weaker than the Fermi interac- 
tion, with Fj «G F = 1.166xl(r 5 GcV~ 2 . If the charac- 
teristic energy scale of these interactions were mj ~ TeV 
then typical couplings are Aj < 1CU 5 for m x ~ 10 keV. 

While x particles may be absolutely stable, extra-weak 
interactions satisfying this bound might allow the decay 
of x particles to SM states on timescales longer than the 
age of the Universe. For instance, for a three-body decay 
we expect a rate like Td ~ Gjm^/(1927r 3 ). A lifetime 
r d ~9x 10 4 (13.7 Gyr) (KT 11 TeV- 2 /G/) 2 (MeV/m x ) 5 , 
orders of magnitude longer than the age of the Universe 
~13.7 Gyr, is possible. We leave the investigation of this 
potentially important signature for future study. 

Minimal Models - There are many possibilities, both 
bosons and fermions, for HHDM phenomenology. As a 
concrete example we have focused here, in analogy with 
SM neutrinos interacting via the Fermi interaction, on a 
model with interactions that decouple when they are well 
approximated by a four-fermion vertex of strength G\. 

A simple model with the right phenomenology consists 
of our dark-matter particle x and another lighter parti- 
cle ip (rriip < m x ) interacting via a spontaneously broken 
U(l) gauge symmetry. If coupled via a massive vector 
boson Ufj_ of mass mjj = m\ with coupling gjj = A, their 
interactions are well approximated by a four-fermion in- 
teraction of strength G\ ~ gfj /rrijj provided m x <C n%u 
and T x f <C mjj. This HS could reheat via direct cou- 
plings of x and ip to the inflaton, or via the hidden Higgs 
sector associated with the broken U(l). A thermal pop- 
ulation of x and ip particles will be produced in the early 
Universe with x accounting for most of the dark matter. 
The ip particles also add to the dark matter density and, 
as they have a longer free-streaming length, might im- 
print detectable features on the matter power spectrum. 

Other alternatives realizing a four-fermion model in- 
clude fermions interacting via a scalar or a vector boson 
with more complex gauge interactions. For instance, the 
interactions of x might be a HS copy of the weak inter- 
actions of SM neutrinos. It would be interesting to find 
well-motivated theories realizing the HHDM scenario. 

Summary - We have shown that a thermalized hidden- 
sector species, with interactions the decouple when rel- 
ativistic, can act as cold rather than hot dark matter 
provided it interacts weakly enough with standard model 
particles to never thermalize with them in the primordial 
plasma. This class of dark matter candidates can account 
for all of the dark matter in the Universe and has an al- 
lowed range of masses between a few keV to a few TeV. 
Signatures of this mass scale may be imprinted on the 
small-scale power spectrum of dark-matter fluctuations 



or, if the dark matter decays to standard model parti- 
cles on cosmological timescales, in the energy spectrum 
of potential decay products. 
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